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Abstract 

The Poisson Nernst-Planck equations for charge concentration and electric po¬ 
tential in a ball is a model of electro-diffusion of ions in the head of a neuronal 
dendritic spine. We study the relaxation and the steady state when an initial 
charge of ions is injected into the ball. The steady state equation is similar to the 
Liouville-Gelfand-Bratu-type equation with the difference that the boundary con¬ 
dition is Neumann, not Dirichlet and there a minus sign in the exponent of the 
exponential term. The entire boundary is impermeable to the ions and the elec¬ 
tric field satisfies the compatibility condition of Poisson’s equation. We construct 
a steady radial solution and find that the potential is maximal in the center and 
decreases toward the boundary. We study the limit of large charge in dimension 1,2 
and 3. For the case of a small absorbing window in the sphere, we find the escape 
rate of an ion from the steady density. 


1 Introduction 

The non-linear system of Poisson-Nernst-Planck (PNP) equations has been widely used 
to study properties of the electric field in local nanodomains such as ionic channels 
[211 [la [ini El E]. It was also used to simulate the equilibration of ions between large 
reservoirs through narrow necks mm- It is also possible to study the effect of interact¬ 
ing ions in ionic channels cite [27119]. 

We use here PNP to study the distribution of charges at the micrometer scale level. In¬ 
deed, the stationary PNP equations with Neumann and no-flux conditions on the bound¬ 
ary of a hnite domain hi, respectively, describe the electrical potential and density of 
charge in hi. They can be reduced to a Liouville-Gelfand-Bratu-type equation for the 
electric potential with however two major differences: hrst, the boundary condition on 
dQ is Neumann but not Dirichlet and second, there is a minus sign in the exponent, 
which is normalized over the domain D. We study here the solution of this equation in 
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spherical symmetry in dimensions < 3 with respect to the (dimensionless) total charge 
A. We construct asymptotic approximations of the solutions for small and large A. The 
one-dimensional case is solved explicitly and it is characterize by a log-singularity at the 
boundary that develop in the large A limit. The explicit solution in two-dimension has 
also a singularity on the boundary. We also obtain a similar asymptotic behavior in 
three-dimension, although the solution cannot be computed explicitly and we provide an 
asymptotic and numerical argument for large A showing again a log-singularity at the 
boundary. We study the voltage change from the center and how it develops a boundary 
layer for large A. The voltage drop from the center to the boundary converges to a hnite 
value as A increases to inhnity. We also apply the analysis of PNP to study idealized den¬ 
dritic spine structure, which consists of a spherical dielectric membrane hlled with ionic 
solution, connected to the dendrite by a cylindrical narrow neck. The paper is organized 
as follow: in the hrst part we study asymptoticall and numericall PNP. In the second part, 
we estimate the current generate in an idealized spine (head connected by a cylinder). 
We derive the current generate by a spine. We show here how the head geometry controls 
the voltage, while the narrow neck radius control the current. 


2 PNP equations in a ball 


We consider the Poisson-Nernst-Planck system in a ball of radius i?, whose dielectric 
boundary dQ is represented as the compatibility condition for Poisson’s equation and its 
impermeability to the passage of ions is represented as a no-flux boundary condition for 
the Nernst-Planck equation. We assume that there are N positive ions of valence z in 
and that there is an initial particle density q{x) in such that 


J q{x) dx = N. 

n 


( 1 ) 


The charge in 17 is 


Q = zeN, 

where e is the electronic charge. The charge density p{x, t) is the solution of the Nernst- 
Planck equation 


D 


ze 


Ap{x,t) + —V {p{x,t)V(j){x,t)) 


dp{x,t) 


D 


dpixA) ze , ^dSixA) 


dt 

= 0 for a; G dfl 


for a? G 17 


p{x, 0) = q{x) for a; G 17, 

where the electric potential in 17 is 4>{x, t) is the solution of the Poisson equation 

A ,/ X zep(xA) , ^ 

A(j){x,t) = -for a? G 17 

EEq 


( 2 ) 

(3) 

(4) 

(5) 
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and the boundary condition 


= -a{x, t) for * e an, (6) 

dn 

where a{x,t) is the surface charge density on the boundary dQ. In the steady state and 
in spherical symmetry 


a{x,t) = 


Q 


47ri?2 


2.1 The steady-state solution 

In the steady state dp/dt = 0 so ([2]) gives the density 


p{x) = N 


exp < 

( ze4>{x)' 

1 

1 kT , 

f 


zecl){x) ) 

> dx 

. kT J 


hence (jS]) gives 


zeN exp 


A(j){x) = 


zecf){x) 

kT 


seq / exp 


ze(j){x) 

kT 


dx 


In spherical symmetry in 


can be written in spherical coordinates as 


(7) 


( 8 ) 


(9) 


0"(r) + 


, , zeN exp ■ 

f ze<p{r )) 

\ 

1 kT J 

'r V ) 
r 

Sd^^Q 1 

[R j 

1 exp < 

ze(j){r ] 

. kT J 

^ r^ ^ dr 


< 0 , 


( 10 ) 


where Sd is the surface area of the unit sphere in The boundary conditions are 

<90(0) _ <90(/?) _ Q 

dr ’ dr SdR‘^~^ 


The inequality in ffTU]) means that 0(r) has a maximum at the origin and decreases toward 
the boundary (see Fig. [3J4)- We can normalize the radius by setting r = Rx for 0 < a; < 1 
and 


u{x) 


ze(l){r) ^ {ze^N 
kT ’ ^ esokT ’ 


( 12 ) 
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to write oa as 


A exp {—M(a;)} 


u"{x) + - — ^u'{x) =- 

SdR‘^~‘^ I exp{—M(a;))} x'^~^dx 

Jo 

m(0)=0, m'(0)=0. 

Incorporating the denominator of the RHS of 0131) into the parameter A by setting 


(13) 


A =/ exp{—M(a:)} ^ dx, 


(14) 


we can write the initial value problem ffT^ as 


u"{x) + - — -u'{x) = — fj, exp {—u{x)} 

X 

m(0) =m'(0) = 0. 


(15) 


First, we show that solutions exist in dimensions 1 < d < 3 only for /i in the range 
0 < /i < yU* for some positive jj*. 


Solution in dimension one 

We solved directly equation [I5] in dimension 1 (see appendix 16.ip and we obtain (see eq. 
1521) that 


u]^{x) = Incos^ 



where I\ is solution of the implicit equation 


h 



(16) 


(17) 


The graph of u]^{x) is shown Fig. [2JA., while the one for versus A is shown in Fig. [2fl. 

We have 0 < /i(A) = ^ and hmA^ooh(-^) = The solution exists u]^ for all A > 0 

and a log-singularity develops at the boundary x = 1 when A —?■ oo. 


Solution in dimension two 

In dimension 2, we obtain the solution in fappendix 16. 2 p 


uTi.x) = log(l 


8/a 


„ 2^2 


X 


(18) 
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where 


h 

/i(A) 

lim yu(A) 

A—^oo 


h 

8 . 


The graph of u\{x) is shown on Fig. [2p, while the one for ^ is on Fig. [2p. u\{x) = 
log(l — develop a log-singularity as A —)■ oo. 


Analysis in dimension three 

The solution of the initial value problem flTT]) in dimension d = 3 can be directly computed. 
We show now that the solution exits for all A, while there is a critical value p*, above 
which, there is no regular solution. Contrary to dimensions one and two, the value of p* 
can only be estimated numerically. We hrst show using a phase-space analysis that the 
solution of equation [TS] is unique when it exists. However it is not possible to use the 
phase-space to study the singularity of the equation. To study the asymptotic explosion 
of the equation, we use an asymptotic argument. Finally, we will study the solution 
numerically. 

Next, we show that the problem flT3|) has a unique regular solution for all A > 0, when 
the solution is hnite. The proof of uniqueness of the solution follows the phase-space 
analysis of flT^ . Indeed, using the change of variables 

s = —logr, u{r) = U{s), v{s) = ^ w = 

ds 

w'{s) = — 2w(s) — U''(s)w(s) = ta(s)[—2 — n(s)], (19) 


which gives 


v'(s) = v(s) — w(s), w'(s) =—w(s)l2 + v(s)], (20) 

and can be written as 

^ ^ -w(2 + v) 
dv V — w 

The phase space of (12(11) contains exactly two critical points. The origin 0 is a saddle point 
and its stable manifold has the tangent T of equation w = 3v. The point Pa = (—2, —2) 
is an unstable node. The initial conditions m(0) = m^( 0) = 0 for the solution of flT^ 
impose lim5^oof^(s) = m(0) = 0 and hms_j.oo fA'(s) = — \\mr^Qru'{r) = n(0) = 0, hence 
the constraints 


lim v{s) = 0, lim w{s) = lim pe = 0. (22) 

S—>-CX) S—>-CX) S—>-CX) 
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Thus the trajectory of the solution of flT^ in the hrst quadrant, which satishes the con¬ 
straints fl2^ . has to be on the separatrix that converges to the saddle point. Choosing 
any value 17(0) gives the value of n(0) = 17'(0) has to be chosen on the separatrix. 

Therefore starting in the hrst quadrant, a trajectory of (|2nil converges to the saddle point 
if and only if it starts on the separatrix with the tangent T. The stable branch at the 
saddle point tends to inhnity as s decreases toward 0. Indeed, the local expansion of (1^ 
near the saddle point is 


m;(u) = 3n + + ..., (23) 

which gives the phase portrait (Fig. [T]). Finally, along the separatrix w'(v) > 0, except 
at the origin, showing that for an initial n(0), there is a unique solution. However, it 
is not possible from the phase-space to study singular solution. Indeed, as we shall see, 
when there is a singularity, because the it occur precisely at the initial value and thus the 
Cauchy problem cannot even start. We conclude the problem has a hnite solution 



Figure 1: Phase-space solution of (|20|) . The separatrix is shown in red, while the other 
trajectories are in blue. 
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and the phase diagram plotted in Fig. [T]ensures that for any initial condition (n(0), ta(0)) 
(when it is hnite) on the separatrix in the hrst quadrant, there is a unique solution to 
(1^ that satishes ([22]). 

A numerical solution of (1T3|) gives the graph Fig. [3E, that is the solution u{x) of (I3|) 
for /i < /i* = 11.2 . The graph in dashed line (/i* = 14) blows up before reaching x = 1, 
while the dash (small point) graph is hnite throughout the interval. To estimate an upper 
bound for fi*, we note that whenever the solution exists for some n near /i*, its asymptotic 
behavior for x close to 1 shows that u''{l) 3> m'( 1) (see the blue graph in Fig. |3|). Indeed, 
to show that under the assumption u''{l) ^ u'{l) the latter inequality is self-consistent, 
we note that near x = 1 the solution of flTdll can be approximated by the solution of the 
simpler problem 


u'\x) = —/iexp {—^(a;)} , 


(24) 


given by 


2 I I 


u{x) ~ log COS“ \ \l 


(25) 


Thus u{x) is hnite in the interval as long as 


71 


< y = 4.934802202 = fi* 


(26) 


and 


u'jx) ^ ^ ^ 

u"{x) ~ y/JjiF 


(27) 


We conclude at this stage that for hxed values of /i below and above, where above the 
latter they blow-up inside the interval 0 < x < 1 (frames A,C,E of Fig. [2|). When fi 
varies with A according to ([HD, the solutions exist for all values of A (frames B,D,F of 
Fig. [2D. Figure EJA-C-E shows that potential drop between the center and the surface of 
the sphere as a function of A for 1 < d < 3. 

In hgure [3l we compare the three dimensional solution obtained numerically with the 
asymptotic expansions in two regimes. We present in appendix 16.31 for A 1, the ex¬ 
pansion u{x) = — A|j^ -|- O(A^) (see eq. [TID. In contrast, for A S> 1, we mention above 
that the approximation u{x) ~ 21n(l — x^), which was relevant near x = 1 can be used 
in the entire domain [0,1]. The analytical approximations (red) are compared with the 
numerical solutions (see appendix [6D. 


The potential differences 

The difference m(0) — u{l) as we shall see in the next section has a physical meaning, as 
it represents the difference of potential between the center and the periphery of a sphere. 
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Figure 2: Numerical solutions u{x) of the initial value problem (fT^ . (A),(C), and (E) 
correspond to different profile values of A in dimensions 1,2, and 3, respectively. The dotted 

curves are solutions that blow-up for x < 1. (B),(D), and (F) are plots of the ratio — in 

h 


dimensions 1,2 and 3, respectively. 








































Figure 3: Asymptotics behavior of the solution u{x). (A) change in the profile 
u{x) for 3 values of the parameter A = 10^, 10^, 10^. (B,C) We present two regimes: for 
A = 0.1 <C 1, we have u{x) = — A|^ + 0(A^) (see eq. UM A 3 > 1 where u{x) ^ 21n(l — 
x'^). The analytical approximations (red) are compared with the numerical solutions (see 
appendix). 


We have in dimension 1, 


I ma(1) -ma(0) I = Incos^ 

where ^ ^ as A —)■ cxo. in dimension 2, 

I ma(1)-ma(0) I = 21og( ^^^g^ ), (29) 

in dimension 3, for A 3> 1, 

I ma(1)-ma(0) I = 21ogln(l -/(A)), (30) 

where the function / is increasing and /(A) —?■ 1 as A —?■ oo. The different curves for 
dimension 1,2 and 3 are shown in fig. 01 In all cases, the large A asymptotic is dominated 
by the log-behavior. 

2.2 Physical implication for the distribution of voltage and charge 
in a dielectric ball 

The distribution of voltage and charge in a dielectric ball can be estimated from the results 
of the previous results by using the dimensional relation [12] in a ball of radius R = 1pm. 
We plotted in Fig. [3J4 the voltage drop for N = 10^, 10^10^ charges. Already for 1000 
charges, there is a difference between the center and the surface of a ball of few milli-Volt. 
This effect could be tested for in the head of dendritic spines. Moreover, the density is 
charge is concentrated at the periphery (Fig. |5j3), leading also to a large field close to the 
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Figure 4: Asymptotics of ma(1) — wa( 0) for dimensions 1,2 and 3. 


boundary (Fig. EP). Consequently most od the charge are accumulated at the boundary, 
as revealed by the plot of the cumulative density of charges 


Q{r) 


N 



f ze(j){r) ) 

dr 

1 kT J 

[K 

/ exp < 

f ze(j){r)' 

1 dr 

1 kT , 


(31) 


We conclude that when the total number of charges is hxed sufficiently high, the charge 
accumulate at the surface. The held is only signihcant close to the surface and thus 
can trap a Brownian charge in such region, while outside a small boundary layer of the 
boundary, the held is almost zero and charge particle experience no drift. This ehect is 
discussed in the section. 


2.3 Scaling laws for the maximum number of charges 

Although we found previously that for a hxed radius, the diherence of potential V (0) — 
1/(1) is bounded as a function of the total number of charge, we shall now show that 
the maximal number of charges increases linearly with the radius of the ball. Indeed, 
introducing the dimensionless radial variable ( = r/R and u\{r) = equation flTUD 

becomes 



Xexp{-Ux/R{C)} 


AttR / exp {-f/A/i?(C)} CdC 


(32) 
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Figure 5: Distribution of (A) the potential, (B) charge and the field (C) and 
cumulative density of charges (D) inside a dielectric ball. 
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with the initial conditions U\/r{Q) = = 0. Now, we solve the initial value problem 


r;(C) + |v';(o = -fiexp{-v;.(0}. v;(o) = v-;(o) = o 

W';(C)=C"exp{-V;(0}, »'„(0) = 0 


and note that 

ux ( t ) = V; (^) . A = 4irpW(l). 


(33) 


(34) 


Thus the number of charges Q in a ball or radius R create the same distribution as a 
charge Q/i? in a ball of radius one, which can be summarized as 


Q{R) = RQiX). 


(35) 


3 Ionic flux in a small absorbing window in a highly 
charged sphere 


We now discuss various consequences of distributing charges close to the boundary, in the 
large charge regime. The first consequence is on the MFPT f{x) from a; G hi, which is 
the solution of the Pontryagin-Andronov-Vitt (PAV) boundary value problem [20] 


D 


ze 


Ar(a;) - —Vt{x) ■ V(j){x) 


= — 1 for a; G n 


dr^x) 

On 


ze , .dchix) 


=0 for X E Qr 


t(x) =0 for X G ria- 


(36) 

(37) 

(38) 


We consider the case of a large field — V0(a3) S> 1 near the boundary |a;| = 1. The profile 
of 4){x) was studied in section 12.11 (see Figures [5]). To study the solution of the PAV 
problem we map the neighborhood of dfla smoothly into the upper half plane 

with coordinates X = {x,y,z), where z = 0 is the image of the boundary, r(X) = f{x), 
and outside a boundary layer near dQa 


V 


d(/)(x) 

On 


= const, 

1 * 1=1 


^{x,y) = 0(a;) 


= const, 

1 * 1=1 


so that Vx,y^{x,y) = 0. The PAV system (I36l) - (l38|) is converted to 

zs 1 

fir„(X) - J^VUX) + K.,f(X) = 


(39) 
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A regular expansion of u{X) for large V gives that to leading order f{X) is a function 
of (x, y) and setting T{x, y) = f (x, y, 0), we find that 

A^^yT{x,y) = (40) 

Thus the MFPT from a; G ^2 to is the sum of the MFPT from x to c?f2 and the 
MFPT form dD. to dD.a on the surface dil. The MFPT to dD. is negligible relative to that 
to dQa- This approximation means that to reach dQa in a highly charged ball a charge is 
hrst transported by the held to the rehecting part dQr of the sphere with overwhelming 
probability and then it hnds dQa by surface diffusion. 


3.1 The current through a small absorbing window in a highly 
charged sphere 

A second conseqnence of the charge distribntions is the control of spine current, inde¬ 
pendently of the voltage. Indeed, the solution T{x,y) of (I40D is the MFPT of Brownian 
motion on a sphere of radius R to an absorbing circle centered at the north-sonth axis 
near the sonth pole, with small radins a = i?sin |. It is given by [25] 


T{x,y) 


“Zf 


log 


sm 


sm 


2 _ 
S ’ 


(41) 


where D is the diffnsion coefficient, 9 is the angle between x and the north pole. Thns 


t(x) = T(x,y). 


(42) 


The MFPT, averaged over the sphere with respect to a nniform distribntion of x is given 
by 

-|- 0(1)^ for (5 1. (43) 

The MFPT for N independent charges is 

fjv = flog - -h 0(1) J for (5 < 1. (44) 

It follows that the cnrrent throngh the small window is given by 

J = ^ = - y —^ for a R. (45) 

2i?2 hog^ + 0(l)j 

We conclnde that once a current enters into a dielectrics ball snch as a spine head, the 
excess of charges Q is hrst pushed toward the boundary and before moving by Brownian 


f = 2R^ (^logi 
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motion to the spine neck. This result shows that the current in a spine head is governed 
by the spine geometry and a key parameter is the radius a of the neck. When there is a 
conservation of charge principle (no leak), the current through the dendritic shaft is the 
same as the one exiting the spine. In that conditions, the spine neck length do not affect 
or modulate the current. 


4 The current in a spine neck under volt age-clamp 
condition 

Determining the voltage drop between the membrane of the spine head and the dendrite 
when a current is flowing from the head to the dendrite remains challenging because the 
cable theory cannot be applied in a system that cannot be approximated by a cable. 
The general scheme for modeling the electro-diffusion in the spine is the PNP model in 
the head and a one-dimensional conduction of ions in the neck. The neck is considered 
a classical ionic conductor. Thus the steady-state PNP equations have to be solved in 
the sphere with boundary conditions implied by the compatibility condition and the flux 
through the neck is determined by the mean hrst passage time (MFPT) of ions from the 
head to the neck, as discussed above. In the case of high charge Q the potential turns out 
to be practically flat throughout the ball with a sharp boundary layer with negative slope 
at the boundary. Thus charge diffuses and is pushed strongly toward the membrane so 
ionic motion is practically conhned to motion on the surface. Due to spherical symmetry, 
the potential is constant on the boundary so ionic motion is free Brownian motion on a 
sphere. At high charge ions interact through the ambient potential that is determined from 
Poisson’s equation in the ball. Therefore they can be assumed independent free Brownian 
particles. The MFPT r of an ion to the small opening of the neck is determined from the 
two-dimensional NET theory (see previous section). Because the flux carried by a single 
ion is q/f, where q is the ionic charge, the number of ions in the spine head is N = Q/q 
and the MFPT tn of any of the N ions is given by 


T 



Thus the current through the neck is 


I = 


Q 


(46) 


and due to charge conservation, it is independent of the length of the neck. If we consider 
the neck to be a parallel-plate capacitor carrying a steady current I, then the voltage 
drop across the neck is simply V = RI, where R is the resistance of the neck, given by 


R 


ksTL^ 
q^nD ’ 


(47) 
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where ks is Boltzmann’s constant, T is absolute temperature, L is the length of the 
conductor, n is the number of ions in the neck, q is the charge of an ion, and D is the 
diffusion coefficient of the solution in the neck [20]. This model is valid as along as the 
voltage is maintained in the spine head. 

5 Discussion, application and conclusion 

We have studied here the solution of the PNP equations in a ball. We estimated the voltage 
in dendritic spines when the voltage in the spine head is maintained. A certain fraction 
of spines receive synaptic connections, essential for neuronal communication. Although 
their functions are still unclear, there are involved in regulating synaptic transmission 
and plasticity [29l[T8],[26l[8| . Interestingly, most of the excitatory connections occurs not 
on the dendrite but rather on spines and the reason is still not clear. The spine shape is 
quite intriguing, made of a head connected to the dendritic shaft by a cylinder. We found 
here that this geometry play a key role: the spine head geometry determines the drop of 
potential, while the current is dehned by the diffusion on the surface and the mean time to 
hnd the entrance of the neck in a two dimensional Brownian motion (see Narrow escape 
time [22|[T3]). In the neck, under a voltage clamp condition, when a constant voltage 
difference between the head and the neck is imposed, the voltage-current relation follow 
a resistance law. Thus the spine geometry dehnes both the capacitance and resistance in 
geometrical terms, a vision that complement previous classical studies [26l[I7l[23|. 

Finally, computing in the transient regime, the change in voltage drop between the 
spine head and the dendritic shaft, requires computing the time dependent PNP equa¬ 
tions. Another open question is to study the influence of the spine head geometry on the 
distribution of charge. Computing the distribution of charges and the associated held in 
non-convex geometry is certainly the most challenging. 


6 Appendix 

In this appendix, we hrst solve analytically the Liouville equation [12] in dimensions one 
and two and in the second part, we describe the numerical methods to compute the 
solution in dimension 3. 


6.1 Solution of the minus sign Liouville-Bratu-Gelfand equation 
in a unitary segment 

Liouville equation in the interval [0 1] is 


- u"{r) 



(48) 
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with initial conditions 


m( 0) = 0 and m^(0) = 0. 


(49) 


This is the classical Cauchy problem. After a direct integration we get with the initial 
conditions 


where 


u’\x) = 


- 1 ), 


(50) 


h = 




(51) 


A second integration gives 


MA((r) = lncos 


A 


Now we self-consistently calculate 


U = / = 


dx 


'0 cos ^\ 


1 / A 

tan \ -. 

A V 2Ja 


Thus Ja > 0 is the solution of the implicit equation 

2 o 


Ja = - tan \ — 


A 


2/a 


(52) 


(53) 


(54) 


The graph of versus A is shown in Figure [2l We have limA^oo ~ specifically, 


Ix 


~ ~ solution fl52|) is shown in Fig. [2] and is regular in the 

entire interval 0 < x < 1 for all values of A. The drop between the extreme points of the 
interval is 


ma(1) - ma(0) = Incos \/^ 


A 


(55) 


and becomes infinite as the total charge increases indehnitely. 


6.2 Liouville equation in dimension 2 

The dimension 2 case can be transformed into the one dimensional case [H] using the 
change of variable 

r = 

u{t) = u{r) — 2t. 
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Equation [13] reduces to 


-uu = Ae-”«)+2> 
where Ix = 2tt e~^^^^rdr and w(t) = u(t) + 2t satisfies 

- Wtt = A—— 


The initial condition are now transform to asymptotic conditions at infinity: 

lim {w{t) — 2t) = 0 


t—>-oo 


lim {w{t) — 2) e* = 0 

t—>-oo 


A first integration gives 


-w{t) 


w'^ e 


2 . 


The solution is 


w(t) = -log(- 


)-2C- 2(, 


(56) 

(67) 

(58) 

(59) 

(60) 

(61) 

(62) 

(63) 

(64) 

(65) 

( 66 ) 

(67) 

A + SvT 

We conclude that ux(r) decreases smoothly and in the limit A —)■ cx), the solution blow-up 
over the entire boundary. 


(Ae2C^+2i - 1)2 
where C is a constant. Finally, we obtain that 

ux(r) = log(l - 

olx 

To close the equation, we shall now compute the integral 


I\= e “^*^’'^27rrdr = 


(1 


A_^2')2 


2'Krdr = 


Svr 


8/a 


8 - A/Ja 


and 


I\ — + Q A 

o 

lim A = 8. 

X—^oo 1^ 

The curve is shown on Fig. |2]and | ma(1) — ma( 0) | in Fig. jH Finally, 


ux{r) = log(l - 


A 


252 


r ) 


A -I- Svr 

ma(1)-ma(0)| = 21og(l - ^ 
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6.3 Regular expansion of solution [15] for small A 

We shall now study the small asymptotic expansion of solution [15] for small A. Using a 
regular expansion, 

u{x) = uq{x) + Mi(a;)A + U 2 {x)X^ + o(A^), (68) 


we obtain using eq. [I5]that uq{x) = 0 and ui is solution of 

1 


Aui 

dui 

dn 


|fi| 


on Vt 




on dVt. 


For i? = 1, 


“‘W = 


with Mi(0) = 0. We conclude that Mi(r) < 0, Thus, 


^0 = -^A + 0(A2). 


The second order term U 2 is solution of 


Ml 




with M 2 ( 0 ) = 0 and = 0- For i? = 1, 

U2{r) = 

Thus, 

uir) =-A 


3r^ 


GdOvr^ 


3r^ 


Svr GdOvr^ 


A2 + 0(A=^). 


(69) 

(70) 

(71) 

(72) 

(73) 

(74) 

(75) 
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